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Abstract 

Stabilization of inter-brane distance is analyzed in 5-dimensional models 
with higher-order scalar kinetic terms. Equations of motion and boundary 
conditions for background and for scalar perturbations are presented. Con- 
ditions sufficient and (with one exception) necessary for stability are derived 
and discussed. It is shown that it is possible to construct stable brane con- 
figurations even without scalar potentials and cosmological constants. As a 
byproduct we identify a large class of non-standard boundary conditions for 
which the Sturm-Liouville operator is hermitian. 



1 Introduction 



Higher dimensional brane models belong to the most interesting recent de- 
velopments in the theory of fundamental interactions. Many models have 
been proposed in which the space-time consists of a 5-dimensional (5D) 
bulk ending at two 4-dimensional (4D) branes. Usually this space-time has 
the structure of a warped product of a maximally symmetric 4D space-time 
and the one dimensional orbifold 5*^/^2 with the branes located at the Z2 
fixed points. The Standard Model fields may propagate only on one of the 
branes called the visible one. Some other fields may live on the second, 
hidden, brane. Of course, the gravity fields can propagate in the whole 5D 
space-time. 

Phenomenological features of such models depend on fields and interac- 
tions other than that of the Standard Model, on the warping, and on the 
distance between the branes. This distance must be fixed in a stable way. 
Such stabilization can not be achieved with only gravity propagating in the 
bulk. A simple mechanism of fixing the inter-brane distance was proposed by 
Goldberger and Wise [1]. The idea is to add a 5D scalar field with some bulk 
and brane potentials. If the background value of that field is not constant 
in the bulk, then the boundary conditions (or in another words: equations 
of motion at the branes) can be fulfilled only if the branes are located at 
appropriate points in the 5th dimension. 

It is not enough to have a background solution with some fixed brane 
positions. It is necessary also that such a configuration is stable against all 
possible small perturbations. From the 4D point of view, the perturbations 
can be describe in terms of Kaluza-Klein (KK) towers of states. The lightest 
scalar KK state is usually called the radion [2]. Tachyonic character of the 
radion indicates instability of a given background. The problem of the radion 
mass, or of the stability of the inter-brane distance, was investigated by many 
authors [3], HJ [5l El [7] . Its relation to inflation was discussed in [8] . Quite 
general criteria for the stability were found in [7]. Generalization of such 
criteria for models with the Gauss-Bonnet interactions was presented in [0]. 

In the present paper we will do the stability analysis for brane models with 
non-standard kinetic terms for the scalar field. Such non-standard kinetic 
terms appear for example in string theory due to the a'- and the loop- 
corrections. Very interesting models with generalized scalar kinetic terms 
were investigated in the cosmological context. Kinetically driven inflation, 
called the k-infiation, was introduced in [TUl [H]. Models of k-essence were 
proposed as another approach to the cosmological constant problem [121 IB] • 
Causality in the context of generalized kinetic terms was discussed by many 
authors (see e.g. [11] and references therein). 
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There is a simple reason to expect that models with non-standard scalar 
kinetic terms may be interesting for the inter-brane distance stabilization. 
Their Lagrangians contain terms with more complicated, than just quadratic, 
dependence on the scalar derivatives. The scalar derivative with respect to 
the 5th coordinate is crucial for the stabilization mechanisms similar to that 
of Goldberger and Wise. This is analogous to the situations in cosmological 
models where the time derivative of the scalar field is crucial. The problem of 
radion stabilization in models with non-standard kinetic terms was addressed 
in [15] but unfortunately the authors used a method which in general is not 
correct and obtained incorrect result^. 

In addition to the bulk (non-standard) kinetic terms we will consider 
also analogous brane-localized ones. Many 5D models with (standard) brane 
kinetic terms for different bulk fields were proposed. Such localized kinetic 
terms were investigated for: pure gravity [ITJ [TBI [IHl [20] , gauge fields [2T], [22l 
ESI El EHl EE], fermions [22l El ES] and scalars ^ EI] 

In section E] we define our model and derive background equations of mo- 
tion and boundary conditions. Analogous equations for the scalar perturba- 
tions are presented in section [SI In subsection 13 . 1 1 we show that the spectrum 
of those perturbations is real. We identify a large class of boundary condi- 
tions for which the Sturm-Liouville eigenvalue problem is self-adjoint. The 
stability conditions are obtained in section [H They are discussed and com- 
pared to that in models with the standard kinetic terms in section [SI Finally, 
section [U] contains our conclusions. 



2 Model and background 

We consider 5D models compactified on the jl^i orbifold with the standard 
gravitational interactions but with non-standard kinetic terms for a scalar 
field $. Two 4D branes are localized at the Z2 fixed points y = yi. The 

^ The authors of [T^ integrated a Lagrangian with fields replaced by their background 
values. They called the result "the effective potential" and looked for the minima of such 
an object. Of course, in general the potential integrated in a given background is not 
equal to the correct effective potential. It happens to be equal in some simple cases, but 
this must be checked case by case by other methods, so the method of integrating the 
potential is practically useless. The authors of [IS] claim e.g. that all models with the 
standard kinetic terms are unstable, what is in clear conflict with the results of many 
previous analyses [3l [3 [71 [HI [3 [16] . 
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action takes the form 



S = jd'xdyV^ !^^R-P{<^,X)-V{<I>) 

-f2~^iy-y^)[Q^H^,X) + U^H^)]\ , (1) 
i=i J 

where 

X = l (V$)^ , (2) 

and 6 is the normahzed Dirac delta satisfying J dy \/—G 6{y — yi) = V— G^*) 
with G^*^ being the determinant of the metric induced on the brane locahzed 
at y = yi (we chose yi < ^2)- The bulk kinetic term is given by some function 
P($, X) depending on the derivatives of $ through the combination X and 
on the scalar field itself. We choose P{^,X) in such a way that it vanishes 
for X = 0. This way V{^) describes the whole scalar contribution to the 
action for constant $. In addition to the bulk interactions, we consider 
brane localized contributions to the scalar kinetic term and to the potential: 
Q»($,X) and respectiveljH . 

The terms in the action ([T]) containing the brane localized kinetic func- 
tions Q*-*-* X) must be treated with special care. Let us discuss now in some 
detail the meaning of an integral containing a product of Q^'^\^,X) and the 
Dirac delta. Writing explicitly the arguments in one of such expressions we 
get 

J d^x J dyv/G(^5(y-y,)g«(^$(x,y),^(G5^(x,y)$''(a;,y) + ...)^, 

(3) 

where prime denotes differentiation with respect to the orbifold coordinate 
y and, in the second argument of Q the ellipsis stand for terms in X with 
derivatives of $ in directions other than y. In brane models the derivatives 
with respect to the orbifold coordinate(s) are usually discontinuous at the 
brane positions. The derivative y), being a Z2 odd function, is exactly 
zero at y = y^. On the other hand, due to the brane sources, the limits 
lim ± y) can be different from zero. The square of the scalar field 

^ There are two kinds of brane kinetic terms considered in the Uterature. Some authors 
assume that such terms involve derivatives with respect to all 5 coordinates (e.g. [21]-[26j) 
while other assume that the derivative in the orbifold direction is not present (e.g. [18j - 
|23| ) . We apply the former approach which seems to be natural when treating thin branes 
as limits of thick ones. Generalization of our results to the case of brane kinetic terms Q 
which do not depend on d^/dy is quite straightforward. 
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derivative, $'^(a;, y), is even under the Z2 symmetry, and can be written 
as a product of sgn^(?/ — y.,) and a smooth function. Usually is 
discontinuous aX y = Ui and, strictly speaking, its integral with the Dirac 
delta localized at yi is not well defined. All expressions of this kind must 
be regularized. Physically, such regularization corresponds to using a thick 
brane and taking the limit of its thickness decreasing to zero. Technically, one 
replaces Siy — y^) and sgn(?/ — y^) with some smooth functions 5e{ii — yi) and 
sgng(y — i/j) satisfying the relation sgn^(y — = 26e{y — yi) and approaching 
the Dirac delta and the signum function, respectively, when e 0. We 
calculate the integrals like ([3]) for regularized expressions and at the end 
remove the regulator taking the limit e — > 0. Thus, we obtain for example 

/ d^a; rdyVG~5{y-y,)^"''{x,y) = [ dV Hm , (4) 

if Va < Hi < Ub- It is not necessary to specify the direction of the limit in (jl]) 
because $'^"(0;, y) is an even function of {y — yi) for any integer n. However, 
the limit itself is necessary because usually |/) is discontinuous at 

y = yi. In the down-stairs approach, one of the limits of integration is equal 
to yi and the r.h.s. of the above equation must be multiplied by 1/2. 

In this work we are interested in warped background solutions with the 
flat 4D foliation described by the ansatz 

ds^ = a^{y){r]i,^dx>'dx'' + dy'^) , (5) 
$ = <P{y). (6) 

The bulk equations of motion for the system described by action ([T]) and 
satisfying ansatz dSHHl) are given by (we use units k = 1) 

{Px<P')' + 3-Px0' - a' {V^ + P$) = , (7) 
a 

!^-2(^:)%W = 0, (8) 

q(^) ^p^^>' + a'{V + P)=0, (9) 



where the subscripts X and $ denote derivatives with respect to the argu- 
ments X and $, respectivel}{§ . 



^ It is straightforward to generalize the equations of motion to the case of any non-flat 
maximaUy symmetric 4D foUation of the 5D background. For example, for the 4D dS 
space-time characterized by the Hubble constant H, the left hand sides of equations ^ 
and ^ should be modified by adding and —6H^, respectively. 
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The boundary conditions for the background can be obtained from the 
full equations of motion resulting from ([T]) with the brane terms taken into 
account. Integrating such equations over an infinitesimal intervals containing 
the brane positions yi one gets 



hm a' = f f/« + / dy6{y - yi)Q^A 



y=yi{y2) 



hm {Fx 00 = ±^ ( + / dy5{y - y.)g« 



(10) 



y=yi{y2) 



where for yi < y2 the upper (lower) signs are to be taken for y = yi {y = 2/2)- 
From now on we use the Dirac delta distribution with the usual normalization 
/ dy6{y) = 1 (in the upstairs approach). 

The above background bulk equations of motion l^lMj and boundary con- 
ditions ( fTOlfTTl) reduce to the know results for the standard kinetic terms after 
substituting P = X and Q^^^ = 0. 



3 Scalar perturbations 

Solving the bulk equations of motion ([319]) and the boundary conditions (ITO] - 
[TTj) one can find possible background configurations characterized by the 
warp factor a{y) and the scalar field (j){y). Not all such background config- 
urations are stable. To check the stability one has to consider all possible 
small perturbations around a given background. Instabilities occur if any of 
the perturbations has a tachyonic character. In this paper we concentrate 
on the scalar perturbations. From the 4D point of view they form an in- 
finite Kaluza-Klein tower of scalars. The state with the lowest (4D) mass 
squared is called the radion. The positivity of its mass squared is a necessary 
condition for the stabilization of the inter-brane distance. 

To find the radion mass we have to investigate the equations of motion 
for the scalar perturbations around the background. Using the generalized 
longitudinal gauge, the scalar perturbations can be written in the following 
way 

ds2 = [(1 + 2Fi) {ri^Ax^dx'') + (1 + 2F2) dy^] , (12) 
$ = + ^3, (13) 

where a and are background solutions depending only on the 5-th coordi- 
nate y, while the perturbations Fj are arbitrary (but small) functions of all 
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coordinates. To find the masses of tlie KK modes of scalar perturbations it 
is enough to consider equations of motion hnear in Fj. 

Contrary to the background equations of motion, for the perturbations 
we obtain non-trivial off-diagonal Einstein equations 



2Fi 



{a'F,)' 



^2 = 0, 



0. 



(14) 
(15) 



They have to be fulfilled in order to stay in the longitudinal gauge. The diag- 
onal Einstein equations, combined with the background equations of motion 
dTHH]), give the third equation for the scalar perturbations: 



□Fi + 4-F; - 4 f-) + - (Px - IxPxx 



-- (Px + 2XPxx) 



b'F. 



b"F. - (h'FL 



0, (16) 



where □ is the 4-dimensional D'Alembertian. The part of the boundary 
conditions linear in the scalar perturbations are quite complicated and reads 

±2 \im [{Px + 2XPxx)F;,] 

+ I 0" [{P<t>x + 2XPm) Fs - {Px + SXPxx + ^X^Pxxx) F^] 



aF. U, 



Tii) 



□ Po 



(17) 



y=yi 



where 6i = 6{y — Ui). The subscript yi at the integrals indicates that the 
range of integration is an infinitesimal interval containing 

The off-diagonal Einstein equations and (fT5l) can be used to express 
two of the perturbations introduced in the ansatz ( |T2lfT3l) in terms of the 
third one. It is convenient to eliminate P2 and P3 and to use the product 
a^Fi as an independent perturbation. We expand it in the 4D modes as 



ikx 



;i8) 



and substitute to eqs. (fT6l) and (fTTl) . Then, the 4D part of the bulk equation 
ffT6l) takes the usual form 



/(m2,fc) + (^^ + f(m2,k) = 0. 



(19) 
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The equation for the "shape" K^2[y) of the KK mode with mass squared 
equal can be written as the Sturm-Liouville equation 



- ipK'^,)' + qK^2 = m'rK^2 , (20) 

where p, q and r are the following functions depending on the background 
3 1 3 



P 



2aP 



x<- 



,2 ' 



1 

Q = - 
a 



2a iPx + 2XPxx] 



(21) 



In the last equality we have introduced a local (y-dependent) "speed of 
sound" 

2 Px 



^ Px + 2XPxx 
The boundary condition 0171) in terms of K^^'^ takes the form 



d 

[hi - CiTn^) ^-^m2 - m^PxK^2 







(22) 



(23) 



where from now on yf stands for y'l or y^. The corresponding limits have to 
be taken for quantities discontinuous on the branes. The 8/ On differentiation 
is in the direction of the outer normal at the boundary, i.e. (— d/d?/) at yi 
and (+d/d?/) at y2- Quantities hi and q are the following functions of the 
background solution and the brane interactions 



aU, 



{i) 



y=yi 



{P<s>x + 2XP$xx) 



T lim {Px + 2XP: 



y^yf 



XX J 



1 

2a 



X 



Vi 



and 



(24) 
(25) 

should be calculated with the same regulariza- 



AU integrals in 
tion as that used in (jlj). 

The square of the radion mass is given by the lowest eigenvalue of the 
equation of motion ( 1201) satisfying the boundary conditions (l23l) . Of course, 
in general it is not possible to find the spectrum of the system fl20l) - fl25l) by 
solving it explicitly. To get some information about the smallest eigenvalue 
we will use methods analogous to those developed for a similar problem in 
[9] (where the corresponding boundary conditions have a form of ( l23i) with 
Cj = 0). But first one has to check whether the differential equation ( l20l) 
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together with the boundary conditions fl2^ constitute a self-adjoint system. 
This is a non trivial problem because conditions fl^^ are unusual and quite 
complicated. The eigenvalue of the equation of motion ( !20l) appears mul- 
tiplying both Krrfl and its normal derivative. In the next subsection we will 
show that our eigenvalue problem is self-adjoint with boundary conditions 
even more general than fl23p . 



3.1 Self- adjoint eigenvalue problem 

Let us consider a differential eigenvalue problem 

Ov = Xv 

for the operator O of the Sturm-Liouville type 

Ov = - [— ipv')' + qv'^ ■ 
The boundary conditions on the interval {yi,y2) have the form 



a?v + afv' + af {Ov) + {Ov)' 



CO 



0, 



(26) 



(27) 



(2^ 



y=yi 



where aj*'' are some constants. The spectrum of our eigenvalue problem is 
real if O is hermitian. In order to prove this one has to find such a scalar 
product (-, ■) for which 

{v,Ou) = {Ov,u) . (29) 

The standard boundary conditions discussed in many mathematical text- 
books have the form of (12H|) with cTg*^ = a^l^ = 0. In such a case, O is 
hermitian in the scalar product (/, g) = J^^ rfg (for simplicity we consider 
real functions / and g). Let us generalize this scalar product by adding some 
boundary term^ 



if, 9) 



y2 



rfg 



yi 



'pff9 + pf{f9)' + pff'9' 



(30) 



with yet unspecified constants . For this scalar product we calculate 

{v, Ou) -{Ov, u) = |p [uv' - vu'] + [v{Ou) - {Ov)u] 

+p? [{v{Ou))' - {{Ov)u)'] [v'{Ou)' - {Ov)'u']] . (31) 

J yi 



^ A simple example of a non-standard scalar product was discussed for example in |27| . 
It was calculated for a canonical kinetic term localized on a branc in a flat background. 
In our notation this corresponds top = 1, q = Q^ r = 1. a\ = 0, (T4 = 0. 
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Introducing three additional constants p^* ^ p^^\ p{\ at each boundary, we 
can rewrite the above equation in the following form 



{v,Ou) - {Ov,u) 



—u 



-u 



pfu - pfu' + pf{Ou) + pf{Ou)' 

'pfv - pfv' + pf{Ov) + pf{Ovy 

'pfu + pfu' + pf{Ou)+pf{Ou)' 
'pfv + pfv' + pf{Ov)+pf{Ovy 



(32) 



where p^2 — PiVi) ~ Pi' ■ operator O is hermitian if the r.h.s. of the 

above equation vanishes for all v and u fulfilling the boundary conditions 
(|28|) . This is the case when each square bracket in fl32l) is proportional the 
square bracket in 



PI 



n\ a 



1 ^1 5 



P2 



Uo (J 



2 ^1 5 



(i) 
-Pl 

pV 



(0 (») 



Pl ■ 

(i) 
P2 



n\ a 



1 ^3 ' 

(i) (i) 



(i) 
P2 

pV 



n a. 



(i) {i 



(33) 



n«a« . (34) 



For generic values of this set of linear equations can be easily solved. At 
each boundary there are 8 equations and 8 independent constants: Pi \ P2\ 
p^^\ p^^\ p^^\ ni \ 77-2'', Pl'' ■ In fact we are interested only in those three, Pi\ 
Ps^ which enter the definition of the scalar product flHUl) . The solution 
reads 



if, 9) 



y2 



rfg 



yi 



P 



/V 



.CO 



y2 



yi 



(35) 



We have shown that the eigenvalue problem flzoj) with the boundary condi- 
tions ( l28i) is self-adjoint. Thus, all its eigenvalues A are real and the eigen- 
functions corresponding to different A are orthogonal in the scalar product 

Let us now use the above result for our k-stabilization mechanism. The 
boundary conditions 0231) have the form of fl28p with 



a 



0, 



^2 



a 



(36) 



with no summation over i. The factors of (—1)* appear because the outer 
normal derivative d/dn was used in (1231) . Substituting ( l36l) into ( l35l) we 
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obtain the following scalar product appropriate to show that the eigenvalue 
problem (I2U1) . (1^ is self-adjoint: 



if, 9) 



-^yi i=1.2 



P'xf9 + Pxc.£{fg) + c^£f£g 



X 



(37) 



The prime at the sum symbol denotes that the boundary contributions should 
be taken only for those boundaries at which Pxbi ^ 0. The reason is that for 
bi = and/or Px = the boundary condition ( |23|) reduces to the standard 
one for which (/, g) = J rfg without any boundary terms (at that boundary). 



4 Stability conditions 

The spectrum of the scalar perturbations in a given background is given by 
the eigenvalues of the Strum-Liouville equation fl^Uj) with the boundary con- 
ditions ( l23l) at the branes. In the previous subsection we have shown that this 
spectrum is real. The most interesting for us is the lowest eigenvalue which 
we identify with the square of the radion mass. The inter-brane distance 
is stable only if this mass squared is positive. In this section we will find 
conditions sufficient for such stability. We will show also when the radion is 
massless and identify some classes of backgrounds which are unstable. 

First we check whether there is a massless mode in the KK tower of the 
scalar perturbations. In such the boundary condition fl23l) at the first 

brane reduces, for nonzero hi, to KQ{yf) = (the case with vanishing bi 
will be considered later). For = 0, the solution of the bulk equation of 
motion ( !20l) . satisfying the boundary condition aX y = yi and normalized to 
-^o(yi) = 1; can be written in quite a simple form 



a\y) 2a\y) 



y 



Ko{y) = ^^7^ - '2/ ^ / dya^(y). (3^ 



Using the background equation of motion ([8]), the derivative of the above 
solution simplifies to 

KM^'-mP4r^yAS)- (39) 



3a{y)a^{yi) 



The boundary condition at the second brane reads 62-^'(z/2^) — 0- The in- 
tegral in eq. (1391) is strictly positive, so this condition is fulfilled only when 
the product b2Px{y2)4>'{y2) vanishes. Repeating the same reasoning starting 
from the second brane, one gets analogous result for the first brane. Putting 
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both cases together, we find that for p and r regular in the bulk, the necessary 
and sufficient condition for existence of a massless mode is 



hhPx{yt)Px{y2Wiyt)<P'iy2) = o. 



(40) 



Conditions sufficient for the stability can be found in the following way. 
Multiplying eq. ( l20l) with K^'^, integrating over the whole 5th dimension, 
and using the boundary conditions ( 123|) we get 



Let us consider first such models for which the background dependent bulk 
functions p, q and r are regular and positive while the brane parameters bi 
are positive and Cj are non-negative. Then, the r.h.s. of fl4ip is positive while 
the l.h.s. is negative for negative and may be divergent for vanishing m^. 
Thus, the condition ( HTl) can be fulfilled only for positive m^. The function 
g = 1/a is always positive. Functions p and r have the same sign as Px and 
[Px + 2XPxx), respectively. They become singular if any of the functions 
Px, {Px + 2XPxx) or 0' vanishes for any value of y. Thus, the following 
conditions 



V,e[,+,,-] r (y) > 0, Px{y) > 0, Px{y) + 2X{y)Pxx{y) > 0, (43) 



are sufficient for the stability of the inter-brane distance (positivity of the 
radion mass squared). By y G [yi,y2] we denote the interior of the bulk, 
yi < y < y2 and the limits y ^ yf and y —> y^. 

Showing that the above conditions are sufficient for stability was quite 
easy. It is much more difficult to check which conditions are necessary. We 
will show now that there must be at least one tachyonic mode if any of 
the functions, (f)', Px or {Px + 2XPxx), vanishes anywhere in the bulk. 
The arguments are similar to those used in [7j and [^. We will compare 
the properties of two solutions of the bulk equation of motion ( 120|) . one for 

= and second for = —M^ in the limit M — > oo. Both solutions 
satisfy the boundary condition at one brane (let us first choose it to be the 
first one located at yi). 



m 



Jy, Y ^ 




h > 0, 



Ci > 



(42) 
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We start with solving the bulk equation of motion fpU|) in the limit of 
large negative = — M^. In the leading order in 1/M, equation (120]) has 
the following approximate solution 



K. 



(y) 



/pCs 



C+ exp 



C+exp 



y 

-Ml Cs ) + exp 

"yi 



y 

-M I c, 
'yi 



-M r cA -C-exp (-M r c, 

Jyi J V Jy\ 



(44) 



(45) 



In the same limit, the boundary condition ( !23l) at the first brane becomes 



y=yi 







(46) 



Because of the M prefactor in fl45|) . for any Ci 7^ and large enough M, the 
above boundary condition can be fulfilled when ~ C~. We choose to 
be positive because later we will compare this solution with Kq normalized 
to 1 at yi. When ci and Pxiyi) have the same sign, the boundary condition 
(HUj) can be fulfilled only when K_M^{yi) and K'_^,p{iii) have the same sign. 
Thus, > C_ and the square bracket in (jHj) does not change its sign in 
the whole bulk. For very large M the first term in (jHj) starts do dominate 
over the second one even for small values oi y — yi (it is slightly bigger even 
at yi) and away from the first brane the solution is approximated by 



K_M2{y) ^ C+^Yy {Px + 2XPxx) exp (^M j 



P 



X 



Px + 2XPxx 



(47) 

Using this solution we can investigate models when some of the conditions 
in (14211431) are not fulfilled. 

It is convenient to define the following function of 



Bo 



m 



{b2 



C2m ) -^K„ 
' an 



m rx 



PxK, 



(48) 



J/=J/2 



It is equal to the l.h.s. of the boundary condition ( l23l) for satisfying the 
bulk equation of motion and the boundary condition at the first brane, and 
normalized to 1 at yi. The spectrum of the KK tower of scalar perturbations 
consists of those values for which B2{m^) = 0. 

Now we check whether the positivity of bi and q are necessary conditions 
for the stability, assuming that all the bulk conditions fj^Hj) are fulfilled. For 
very lage negative the boundary function B2 at the second brane is dom- 
inated by the term proportional to K'_j^,^2- From eq. ( 1471) and the discussion 
before it, it follows that 



sgn [B2{-M^)\ = sgn [M'^c^K^.^.^y^)] = sgn [ca] 



(49) 
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On the other hand, for the solution Kq given by flHKI) and flH^ we get 

sgn [52(0)] = sgn [62^^0(2/2")] = , (50) 

where we used the fact that Kq is always positive when the inequalities 
are fulfilled. Comparing fHUj) with fl3U|) . we conclude that there must be at 
least one negative eigenvalue when the parameters 62 and C2 have opposite 
signs. For 62C2 < 0, the function B2{'w?) has different sign for m? = and 
for large (enough) negative m^. There must be some negative for which 
B2 vanishes because the solutions of (1201) change continuously with m^. 

Repeating the above reasoning but starting from the brane at y2, we 
obtain an analogous condition for parameters bi and ci. Thus, the conditions 

6ici>0, fc2C2>0, (51) 

are necessary for the stability. 

Now we investigate the stability conditions for the bulk quantities ( l43l) . 
The solution (H7j) for large negative vanishes at a point at which or 
Px vanishes. It must change sign there because from fl^Ul) it follows that 
K and K' can vanish at the same point only for trivial solution vanishing 
everywhere. Thus, for very large negative the function K{y) vanishes 
close to the point where Px4>' is zero. On the other hand, from ( l38l) and ( l39l) 
it follows that Kq is positive for all y. So, there must be some negative 
for which Kf^2 has a zero point but is nowhere negative. It is easy to see 
that such a zero point must be at the second brane, y = y2, and that the 
derivative of K^2{y2) is negative. In such a situation 

sgn [B2{m^)\ = sgn [(62 - C2m'^)K'~2{y2)] = -sgn[62] , (52) 

where the last equality follows from the condition 0511) . Comparing fl50l) and 
(l52l) we find that there must be some negative mode with the eigenvalue 
satisfying < < for which B2{m^) = 0. The radion is tachyonic if (p' 
or Px vanishes in the bulk. 

The above arguments are rather complicated but the result is quite in- 
tuitive. We consider backgrounds for which Px4'' vanishes at some yo < 1/2 
in the bulk. For any such background Ko{y) defined in ( !38|) is a zero mode 
in a model restricted to the interval [yi,yo\. It is quite natural that the KK 
states becomes lighter when the compact space becomes bigger. So, with 
a massless mode on [?/i,yo] there should be a tachyonic one on the bigger 
orbifold [yi,y2]- 

Equation ( H71) can be used to show that also {Px + 2XPxx) should be 
strictly positive. If it is not, there are two possibilities depending on how 
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fast it approaches zero. If the integral in fjTTl) is finite then K_m2 vanishes 
because of [Px + '^XPxx) in the prefactor and a reasoning similar to that 
for the case of vanishing Px<P' niay be applied to prove the existence of 
at least one tachyonic mode. On the other hand, a divergent integral in 
(H71) indicates the breakdown of the perturbativity assumption. This is not 
surprising. Vanishing (Px + 2XPxx) corresponds to infinite speed of sound 
while negative (Px + 2XPxx) gives negative square of the speed of sound 
(for positive Px, which is anyway necessary for the stability). In both cases 
one should expect strong instabilities. 

We showed above that the conditions ( l43l) on the bulk quantities are not 
only sufficient but also necessary for the stability. We were not able to prove 
the same for the brane conditions fjl21) . If one of them is fulfilled then the 
other has also to be fulfilled. The only possible loophole occurs when both 
conditions ( l42i) are violate, namely when 61 < and ci < or when 62 < 
and C2 < 0. However, these are not very appealing possibilities. Parameters 
Ci are proportional to the integrals ^iQx ^'^^ can be negative only for 
localized brane kinetic terms very different from the standard one. 

5 Discussion 

With the results presented in the two previous sections we can investigate 
how the stabilization of branes is influenced by the presence of non-trivial 
scalar kinetic terms in the bulk and/or on the branes. Such terms change 
the background configurations and the spectrum of the scalar perturbations. 
We start the discussion with the background. 

Combining eqs. (IHl) and (l9l), the dynamical equation describing the change 
of the warp factor can be written as 

a" 

3— + a2(y + P) = 0. (53) 
a 

The source for the change of the warp factor a{y) is the full "matter" La- 
grangian density {V + P) irrespective of whether the kinetic part is standard 
or not. The modification of the scalar equation of motion given in ([7]) 

(Px0')' + 3- (Px0') - (Ke, + p.,) = , 
a 

is twofold. First, similarly as in the case of the warp factor, the role of the 
potential in this equation is played by the full non-gravitational Lagrangian 
density. Second, it seems that a natural variable to describe the change of 
the scalar background is the product Px(t>' and not (p' itself. The equation of 
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motion for this generalized variable Px'P' looks formally the same as that in 
the standard theory (derivative of the full Lagrangian as a source and 3a'/ a 
as "friction"). Thus, as compared to the standard theory, for the same local 
non-gravitational energy density and the warp factor slope, the scalar field 
changes faster (slower) if Px is smaller (bigger) than 1. Of course this is only 
a qualitative feature and in most of the cases any quantitative corrections 
can be found only by numerical calculations. 

The positions of the branes are determined by the boundary conditions. 
The modifications to the boundary conditions (ITOl) and (fTTl) are analogous 
to those in the bulk background equations. Namely, not only the potentials 
but the full Lagrangians localized at the branes determine the jumps of a' 
while their derivatives with respect to $ determine the jumps of Px'P' ■ 

Usually in Randall-Sundrum type models, the warp factor changes mono- 
tonically in the bulk, so its derivative has the same sign for all y. Thus, be- 
cause of opposite overall signs in the boundary conditions ( ITOi) at two branes, 
one brane must have positive tension while the second one must have nega- 
tive tension. To check the signs of the brane tensions in the class of models 
considered in this work we rewrite eq. ([H]) in the following form 

-)i = -^Px<P'' . (54) 

In the previous section we showed that the stability of the model requires 
that Px'P' is everywhere non-zero. Thus, the r.h.s. of the above equation is 
always negative. The ratio a'/a^ always decreases and the warp factor aiy) 
can not have a minimum in the bulk. Because of that, it is not possible to 
construct a stable model with two positive tension branes. At least one brane 
must have a negative tension: 



mm 



ij/i 



U^\+ / ^^Q^''A <0. (55) 



In all stable models 0(?/) must be a monotonic function (0' can not vanish) 
and Px can not change sign. The limit of the product Px'P' has the same 
sign at both branes. Thus, it follows from the boundary condition fill I) that 

(u'i^ + / 5igi')) ■ ff/f + / 52Qt^\ <0. (56) 

We turn now to the stability conditions. One of them is the positivity of 
hi parameters defined in fl23|) . Using the background equation of motion ([7]) 
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the last term in the definition of bi fl^^ can be rewritten as 



T hm^ {Px + 2XPxx) 



( 



4>' 



o! 



a 




(57) 



where we used the outer normal derivative introduced in eq. fj23l) . The first 
term in the last square bracket of the above equation gives negative (positive) 
contribution to the b parameter on the positive (negative) tension brane. 
So, positivity of b at the positive tension brane is more difficult to achieve. 
Stability is improved when, close to the brane(s), P^x has the opposite sign 
and (V$ + P$) has the same sign as the normal derivative of the scalar field 
d(j)/dn. There is another term in the definition of b which depends on the 
bulk background, namely / 0"(P$x + 2XP^xx)- Its sign depends on the 
background and on the details of the generalized bulk kinetic function P. 
Non-trivial ^-dependence of P can be, at least in some cases, used to increase 
the radion mass. Finally, large enough values of the second derivatives of the 
brane kinetic terms Q*^*-* may be used to make 6j positive. 

The second stability condition in 0421) can be quite easily fulfilled. For 
example: q given by eq. (l25l) vanishes if there is no kinetic term localized 
on the i-th brane and it is positive when such localized term is not much 
different from the standard one Q*^*-* = X. 

Models with non-standard bulk and/or brane scalar kinetic terms are 
quite complicated and usually only performing numerical calculations one 
can find the background fields and check their stability against small pertur- 
bations. Nevertheless, it seems viable that stable solutions can exist also in 
models without any scalar potentials or cosmological constants. The kinetic 
terms alone may have structure rich enough for configurations with stabilized 
inter-brane distance. This is similar to the situation in models proposed in 
pUl [TTj in which infiation was realized without any scalar potential. 

Let us discuss what properties the generalized kinetic terms should have in 
order to support stable brane configurations. Conditions on the bulk kinetic 
function P are rather weak. It is enough that eq. ([9]) can be fulfilled for some 
Uq and positive values of 0'^, Px and {Px + 2XPxx)- Then, the dynamical 
equations ([7]) and ([H]) can be used to extend the solution to y ^ y^. The bulk 
stability conditions (H3l) are fulfilled at y^, so they are fulfilled also in some 
finite interval in the 5th coordinate. Any two points in this interval may be 



17 



used to locate the branes. Of course, this is possible only when the brane 
kinetic terms have appropriate properties. 

Restrictions on the brane kinetic functions Q^*^ are quite strong if we 
want the branes to be stabilized at given positions in a given background. 
First of all, from eq. ( l56l) it is obvious that without brane potentials it is 
necessary that Q*-*-* have some non-trivial ^-dependence. In addition, it 
follows from that at least at one of the branes the kinetic term must 
give a negative contribution to its tension. This does not a priori mean 
that the system becomes unstable. Of course, we want the energy to be 
bounded from below, so the kinetic term at the second brane Q^'^\^^X) 
(we call "second" that brane at which the expression in fl55|) is minimized) 
should give negative value of J S2Q^^^ only for some range of values of its 
argument^. More specifically, each Q*-*-* must satisfy two equalities ffTUj) and 
(ITT]) and two inequalities (H^ . The values of / 6iQ^^'' and / SiQ^^^ necessary 
to fulfill the background boundary conditions depend on the details of a given 
background but their signs are determined by (l55l) and ( l56l) . All boundary 
and stability conditions on the brane kinetic functions may be written in the 
following form 



ii) 



6 da 
a? dn 



5^Q 



a dn 



(58) 
(59) 



/ 5^Q5^ >0, (60) 

S..Qg>-^,, (61) 

where bi is the r.h.s. of ( 12^ with Q*-*^ set to zero. It is possible to fulfill all 
the above conditions for example with the brane kinetic functions of the form 

Q W = ^« + L» ($)X2 . (62) 



The most difficult part is to satisfy simultaneously conditions ( 1581) and ( 1601) 
at the second (negative tension) brane. Using eqs. (l58l) . ( 1601) and (jll), one 



^ One should remember that in general Q^*) and J SiQ^"^^ are not just proportional to 
each other and can have negative values for (slightly) different regions of the parameter 
space. This is caused by the regularization procedure discussed in section [2] 
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can find a lower bounds on L^*) 



' yi a'^ on 



which can be translated to an upper bound on i^^'^ 



(63) 



+ 



18 9a 
2 dn 



a 



< 



Vi 



180 da 



dn 



(64) 



Vi 



At the positive tension brane L^^^ (0(7/2)) niust be positive and big enough 
while ii'^^'' (0(7/2)) must be negative (with value related to the value of L*^^)). 
Thus, without scalar potentials it is not possible to construct a stable model 
with a positive tension brane if the corresponding K^'^'> is always positive. 

Some ^-dependence of i^^*^ and/or L*^*) is necessary to fulfill conditions 
(l59l) and (1611) . The background boundary condition (l59l) takes the following 
form 

2P.d, ^^^^ 



(') _L T "^2 



9n 



Vi 



In all stable configurations, the r.h.s. of this equation has opposite signs on 
the two branes (because 0' can not change sign). So, there are no consistent 
brane models without potentials if all first derivatives of K'^'^'^ and L^*) have 
the same sign. 

The stability conditions (161 p for the brane kinetic functions (1621) may be 
rewritten as 



lim 

y^vt 



g-f'-^Cli^ ^^^^ 



<i>X 



P 



X 



3 * 



(0 



1 



X' 



APx \ -K^'^X + —L^'^X^ + - 



\Kt^X + \L 



X2 



> 



(66) 



Vi 



Some of the terms on the l.h.s. of the the above expression may be negative 



but they can be compensated by large enough value of Kl^X/?, + L^^XVS. 

It is clear that it is possible to choose functions i^^*^ and L*^*^ which satisfy 
all the above boundary and stability conditions for a given background. So, 
models in which the inter-brane distance is fixed in a stable way can be 
constructed even without any scalar potentials or cosmological constants. 
The brane induced kinetic terms may have a relatively simple form Q = 
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KX + LX^ if the functions K and L are generic enough|§. It would be 
interesting to check whether any higher order kinetic terms predicted for 
example by string theories have an appropriate structure. 

6 Conclusions 

We considered 5D brane models with bulk and brane scalar kinetic terms 
generalized to some functions oi X = (V0)^/2 and the scalar field itself. 
The background equations of motion and boundary conditions have struc- 
ture similar to the case with standard kinetic terms. There are two kinds 
of modifications. First: the scalar potential is replaced by the sum of the 
potential and the kinetic term. Second: derivatives of the scalar field are 
multiplied by derivatives of the bulk kinetic term with respect to X. 

Stability of background configurations has been checked by analyzing 
the spectrum of small scalar perturbations. A given background with fixed 
branes positions is stable only when all the masses squared in the spectrum 
are positive. The bulk equation of motion determining the shape of the 
KK modes of such perturbations was written in the Sturm-Liouville form. 
The corresponding boundary conditions have rather complicated form. They 
may be expressed in terms of four parameters (two for each brane), hi and 
Cj, determined by the background and by the bulk interactions described ef- 
fectively by some potentials and generalized kinetic terms. The boundary 
conditions depend also on the eigenvalues and this dependence is more com- 
plicated than in models with standard kinetic terms. We have shown that 
our eigenvalue problem is self-adjoint with those complicated boundary con- 
ditions. We identified even larger class of boundary conditions for which the 
Sturm-Liouville operator is hermitian. 

The eigenvalue-dependence of the boundary conditions makes the sta- 
bility considerations more difficult. Sufficient conditions for the stability 
are: 6j > 0, q > at each brane and the positivity of bulk functions Px-, 
{Px + 2XPxx) and 0'^ for all values of the 5th coordinate y. If Cj > then 
the remaining conditions are not only sufficient but also the necessary ones. 
This changes when any of the q parameters is negative. It seems that it may 
be possible to have stable configurations with negative both bi and ci (or 62 
and C2). The lowest KK mode, the radion, becomes tachyonic when any of 
the quantities 6jCj < or any of the quantities 0'^, Px or [Px + 2XPxx) is 
not strictly positive. 

^Of course, one fine tuning of parameters is necessary as in all models with flat 4D 
foliation. 
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We have shown that stable brane models may be constructed without bulk 
and/or brane potentials and cosmological constants. This may be achieved 
for example when the brane localized kinetic terms take the form Q^*^ = 
K(^(^)X + Conditions for the functions K^'\^) and 

have been found. 
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